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The explicit expressions for the non-singlet DIS structure functions Fi and < 71 , obtained at small 
X by resumming the leading logarithmic contributions to all orders, are discussed and compared in 
detail with the DGLAP evolution for different values of x and . The role played by the DGLAP 
inputs for the initial parton densities on the small-x behavior of the non-singlet structure functions 
is discussed. It is shown that the singular factors included into the fits ensure the Regge asymptotics 
of the non-singlet structure functions and mimic the impact of the total resummation of the leading 
logarithms found explicitly in our approach. Explicit expressions are presented which implement 
the NLO DGLAP contributions with our small -x results. 
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I. INTRODUCTION 

The non-singlet components of the structure function Fi (= and of the spin structure function gx (= have 
been investigated in great detail in deep inelastic scattering (DIS) experiments. The standard theoretical framework 
is provided by DGLAP0]. In this approach, any of f^^{x,Q^),g^^{x,Q^) can be represented as a convolution of 
the coefficient functions and the evolved quark distributions calculated in LOQ, NLOgj and NNLO[3. Combining 
these results with appropriate fits for the initial quark distributions, provides a good agreement with the available 
experimental data (see e.q.Q and the recent paper 0) 

On the other hand, the DGLAP evolution equations were originally applied in a range of rather large x values, 
where higher-loop contributions to the coefficient functions and the anomalous dimensions are small. According to 
the results of Refs HQ, such corrections are becoming more and more essential when x is decreasing and DGLAP 
should not work so well at a; I. Nevertheless, it turns out that DGLAP predictions are in a good agreement with 
available experimental data, leading to the conclusion that the impact of the higher-order corrections is negligibly 
small for the available values of x and those corrections may be relevant at asymptotically small x only. 

In the present paper we use our previous results Hi to show that the impact of the high-order corrections on the 
and x -evolutions of the non-singlet structure functions is indeed quite sizable at available values of x. We also 
show that the reason for the success of phenomenological analysis of the data based on DGLAP at x < 10“^ is related 
to the sharp x -dependence assumed for the initial parton densities, which is able to mimic the role of high-order 
corrections. 

The paper is organized as follows: In Sect. 2 we introduce our notations and discuss the difference of our approach 
with DGLAP. In Sect. 3 the DGLAP and our predictions for the small-x asymptotics are explicitly compared. It 
allows us to clarify the role played by the singular terms in x in the DGLAP inputs for the initial parton densities. In 
Sect. 4 we give a detailed numerical comparison of our results with DGLAP, and discuss also the role played by the 
regular part of the initial parton densities. As our approach deals with small x only, in Sect. 5 we suggest a possible 
method to combine DGLAP with our approach in order to obtain equally correct expressions for both large and small 
values of x. Finally, Sect. 5 contains our conclusions. 
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II. COMPARISON BETWEEN DGLAP AND OUR APPROACH 

As the DGLAP -expressions for the non-singlet structure functions are well-known, we discuss them briefly only. 
In this approach, both and gi^{x,Q'^) (= fDGLAp{x,Q^ j) can be represented as a convolution 

fDGLAp{x,Q‘^) = [ —C{x/y)Aq{y,Q'^) (1) 

Jx y 

of the coefficient functions C{x) and the evolved quark distributions Aq{x, Q^). Similarly, dAq{x, Q"^)/ can be 

expressed through the convolution of the splitting functions and the initial quark densities (5q(a; « 1 , ~ /i^) where 

is the starting point of the -evolution. It is convenient to represent fDGLApix,Q'^) in the integral form, using 
the Mellin transform: 


foGLApix, 


di ) 

— {l/x)‘^C{uj)6q{u})exp 




( 2 ) 


where C{uj) are the non-singlet coefficient functions, 7 (w) the non-singlet anomalous dimensions and 6q{uj) the Mellin 
transforms of the initial non-singlet quark densities. In the x- space, the standard DGLAP inputs for 5q{x) for the 
non-singlet parton densities (see e.g. Refs. ii) include terms which are singular when x ^ 0 , as well as regular 
ones for all x. For example, the ht A of Ref-jJ is chosen as follows: 


Sq{x) = Nyx “(I — x)^(I-I-= Nyx °‘(j){x) , (3) 

(j){x) = (1 - x)^{l + jx^) , 

with N, r] being normalization factors, a = 0.576, (3 = 2.67, 7 = 34.36 and 5 = 0.75. As the first term NrjX~°‘ in the 
rhs of Eq. © is singular when x —*■ 0 whereas the second one, (/'(x) is regular, we will address them as the singular 
and regular parts of the fit respectively. Obviously, in the lo -space Eq. © is a sum of pole contributions: 


k^l 


5q{uj) = Nr] (a; — a) > rak ((u) k — a) ^ k 1 — a) 




( 4 ) 


with ruk = f3{P — l)..(/3 — k + l)/fc!, so that the first term in Eq. Q) (the leading pole) corresponds to the singular 
term x““ of Eq. © and the second term, i.e. the sum of the poles, corresponds to the interference between the 
singular and regular terms. In contrast to the leading pole position co = a, all other poles in Eq. 0 occur at negative 
values, because fc — a > 0. As usual, when co is integer: to = n = 1,2,.., the integrand of Eq. © gives the n-th 
moment of /, (= Tn{Q^))- The standard DGLAP -technology of calculating f(x,Q‘^) corresponds to first calculating 
TniQ^) and then reconstructing f{x,Q^) from the moments, choosing appropriate forms for Aq. Presently C{uj) for 
the non-singlet structure functions are known with two-loop and three- loop [3 accuracy. 

In order to make the all-order resummation of the double-logarithmic contributions to /^'®(x, Q^) and ^, 
i.e. the contributions ~ {as/io'^Y fo 7 ( 0 ;) in Eq. |(2I), an alternative approach was used in Refs, p, by 

introducing and solving infrared evolution equations with fixed Og. This approach was improved in Refs. 0, where 
single-logarithmic contributions were also accounted for and the QGD coupling was running in the Feynman graphs 
contributing to the non-singlet structure functions. In contrast to the DGLAP parametrization Og = ag(fcp, we used 
in Refs. 0 another parametrization where the argument of as in the quark ladders is given by the time-like virtualities 
of the intermediate gluons. Arguments in favor of such a parametrization were given in Ref. Q- In particular, it was 
shown that it coincides with the DGLAP -parametrization when x ~ 1, but it differs sensibly when x <C I. Instead 
of Eq. ©I Refs. 0 suggest the following formulae for the non-singlet structure functions 





^{l/xYC^^l{u})5q{uj) exp (iJ^j(a;)?/) 


(5) 


Integration over lo in Eq. 13 was performed in Refs. however with simplified assumptions for ols 
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rJ( 

— (w)( 5 ( 7 (w) exp 


( 6 ) 


with y = ln((3^//i^) so that y?' is the starting point of the -evolution. The new coefficient functions are 

expressed in terms of new anomalous dimensions : 


r<(±) _ 




(7) 


and the new anomalous dimensions H^g (oj) which account for the resummation of the double- and single- logarithmic 
contributions are 


= (1/2) L - v'a;2-il±(a;) 


( 8 ) 


where 

= (47rCr(l + u;/2)A(u;) + D<^^^u;))/(27r^) . (9) 

Finally and A(u>) in Eq. 0 are expressed in terms of p = ln(l/x), b = (33 — 2n/)/127r and the color factors 

Cf = 4/3, TV = 3: 


D^^\uj) = 


2Cf 


ln(^^) 


77 


P + V 


T 


-{p + pY + p- 


( 10 ) 


1 r 77 r°° dpe 

51.772-1-772 Jq (p-I- 77)2-I-7r2. 


( 11 ) 


A{uj) is the Mellin transform of as with the time-like argument as{k'^) = l/(61n(—/c^/A^)). The comparison of 
Eq. 0 to Eqs. 00 clearly shows the difference between DGLAP and our approach: 

First, our approach accounts for the double- and single- logarithmic contributions to all orders in the QCD coupling 
whereas in the DGLAP -expressions the coefficient functions and the anomalous dimensions are obtained up to 
fixed-order accuracy (two and three loops). 

Then a second difference is related to the treatment of as- In DGLAP, the argument of ag in any Feynman graph 
is k'^, with k being the momentum of the ladder quarks, and therefore as is not involved in the Mellin transform. 
This leads to the dependence of the non-singlet structure functions with the exponent ~ InlnQ^. In our approach 
as depends on time-like virtualities of the ladder gluons and therefore, it is involved in the Mellin transform. 

The inclusion of as in the Mellin transform explicitly shows that the factorization between the longitudinal and 
transverse spaces used in DGLAP is not valid for small x. However, as shown in Ref. |^, our approach and DGLAP 
converge when x ^ 1 where the factorization does hold. The parametrization of as{k‘^) in the form of Eq. has 
recently been used in Refs. m for calculating the double- and single- logarithmic corrections to the Bjorken sum 
rules. 


III. THE ROLE OF THE SINGULAR FACTOR FOR Sq IN DGLAP FITS 

In the first place let us study the small-x asymptotics of the non-singlet structure functions. When x —> 0, one can 
use the saddle point method in order to estimate the integrals in Eqs. 15I6II and derive much simpler expressions for 
the non-singlet structure functions. When the initial parton densities Sq{x) are regular in x (e.g. given by (/(x) only 
in Eq. 0 ), they do not contribute to the saddle point position, which is: 


with 


= /^[l + (I-(7(±))2(77/2 + I/\/(^)V(21n2^) 


= dVB(±)/du 


( 12 ) 

( 13 ) 




















and y = ^ 

It immediately leads to the Regge asymptotics for the non-singlets: 
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( 14 ) 


with 




■1 1/2 , 


2(l-g(±))VR(±) {y/2 + l/VW^) 


7ri/2 In^/^ ^ 

For Q2 ^ 150 GeV^, y < 2 //bW, so ojg^^ do not depend on y. Therefore Eq. d can be rewritten as: 


(15) 


^ e2<5g(4+))c(+)T(+)(C), g^^ix,Q^) ^ el5q{uj^g-^)c^-^T^-\0 


(16) 


with 


and the intercepts = 0.38 and ojq 


(-) 


T(±)(^) = ^“o*’/2/1^3/2^ 

= 0.43 (see 0 for details). The factors 

r -n 1/2 

c(±) = 2(l-g)/(7rVR(±)) 


(17) 


(18) 


do not depend on y. Eq. d predicts the asymptotic scaling for the non-singlet structure functions: in the region 
<C 150 GeV^, with T^^'> depending on one argument ^ instead of x and independently. In the opposite case 
when > 150 GeV^, the scaling is valid for /y and /y- The small-x dependence of Eq. HI 611 is confirmed 
by various analyses of the experimental data in Refs. la, while the Q^- dependence of these formulae has not been 
checked yet. 

On the other hand, under the same hypothesis that 6q{x) is regular (e.g. given again by (j){x) in Eq. Q), the 
small-x asymptotics of /dglap is well-known (we drop the inessential pre-factor here): 


foGLAP ~ exp y {{2CF)/TTb) ln(l/x)((ln(Q2/A2))/(ln(ii2/A2))^ . (19) 

Obviously, this behavior is different and less steep than the Regge asymptotics of Eq. d. However when the singular 
part in 5q is also included in the DGLAP input of Eq. ®, the DGLAP asymptotics for the non-singlet structure 
functions is controlled by the leading singularity u} = a = 0.576 of the fit in the integrand in Eq. 0 so that the 
asymptotics of foGLApix, Q^) is also of the Regge type: 

JSSlap - (e2/2)C(a)(l/x)“((ln(QVA2))/(ln(^VA^)))''^“^^'’ , (20) 

with b = (33 — 2nf)/12TT. The comparison of Eq. (I16II and Eq. (I20II shows that both DGLAP and our approach 
lead to the Regge asymptotic behavior in x for the non-singlet structure functions. In particular, Eq. 0 makes the 
DGLAP prediction be more singular in x than ours, Also they differ as far as the predictions for the -behavior are 
concerned. However, it is important to stress that our intercepts uj^ are naturally obtained by the total resummation 
of the leading logarithmic contributions, without assuming any singular input for Sq, whereas the DGLAP intercept 
a in Eq. d is generated by the phenomenological factor x of Eq. which makes the structure functions grow 
when X decreases, and mimics in fact the total resummation^. In other words, the role of the higher-loop radiative 
corrections on the small-x behavior of the non-singlets is, actually, incorporated into DGLAP phenomenologically, 
through the initial parton densities fits. It also means that the singular factors can be dropped from such fits when 
the coefficient functions account for the total resummation of the leading logarithms and therefore regular functions 
of X only can be chosen for the initial densities Sq. 


2 We remind that our estimates for the intercepts were confirmed by several analyses!^ of the experimental data 
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IV. THE ROLE OF THE REGULAR PART IN THE INITIAL PARTON DENSITIES 


As discussed in Sect. II, the regular part in the initial quark distribution in the DGLAP approach, e.g. (j){x) = 
N{1 — x)^(l + 'yx^) in Eq. ®, corresponds to the sum of pole contributions (t>{uj): 


OO 

= N'^mt + A;)“^ + 7 ( 0 ; + k + 1 )“^^ 

k=l 


( 21 ) 


We have introduced here the normalization factor N. Obviously, (jj —f N when x —> 0. An input Sq{x) = N means that 
the shape of the initial quark distribution does not depend on x. In order to study the influence of the initial parton 
densities in more detail, we compare numerically our results for two situations: (a) a regular fit 6q{x) = (j){x) is 
used, (b) 5q{x) = N. We define Ri{x) as follows: 


Ri{x) =g^^[Sqix) = (l)ix)]/g^^[Sqix) = A] (22) 

at = 10 GeV^. For all numerical studies in the present section we choose the starting point of the Q -evolution 
p. = 1.5 GeV and assume A = 0.15 GeV and n/ = 4. The results for i?i are plotted in Fig. 1 (curve 1). It shows that 
one can approximate 4>{x) by a constant distribution for x < 0.05 whereas at larger values of x the x -dependence of 
the initial quark distributions is essential. 

The comparison of our results @1 with DGLAP when the regular fit (j>{x) is used is also shown in Fig. 1 (curve 2) 
for the ratio R 2 defined as 


R 2 ix) = g^^ [Sq{x) = (j3{x)] /gi DGLAP = (jj^x)] (23) 

where we choose again = 10 GeV^. This shows that when the input for <5q does not have a singular term, the 
difference between our results and DGLAP grows fast for x < 0.05. On the contrary, when the singular term is 
included into the DGLAP input, the DGLAP results grow faster than ours. This is shown in Fig. 1 (curve 3) for the 
ratio R 3 : 


Rsix) = gi^ [d< 7 (x) = <()(x)] /gi DGLAP = x “<()(x)] , (24) 

with = 10 GeV^. A fast growth of foGLApi^^i^) = x~°'(j){x)] for small x also contradicts the analyses of 
experimental data obtained in Refs. |15| . 

Our results 15161) differ slightly from DGLAP also for the -evolution, depending on x. This is shown in Fig. 2 
where the ratio R{Q^): 


RiQ"^) = 9i^ [Sq{x) = <()(x)] /gi DGLAP [^qix) = 4>(x)] (25) 

is plotted for x = 0.01 (curve 1), x = 0.001 (curve 2) and for x = 0.0001 (curve 3). Fig. 2 also suggests that the 
difference between our and DGLAP predictions for the -evolution slowly grows with increasing Q^. 

V. IMPLEMENTATION OF DGLAP WITH OUR HIGHER-LOOP GONTRIBUTIONS 

We have shown that Eqs. account for the resummation of the double- and single logarithmic contributions to 
the non-singlet anomalous dimensions and the coefficient functions that are leading when x is small. However, the 
method we have used does not allow us to account for all other contributions which can be neglected for x small but 
become quite important when x is not far from 1. On the other hand, such contributions are naturally included in 
DGLAP, where the non-singlet coefficient function Cdglap and anomalous dimension 'Jdglap are known with the 
two-loop accuracy^: 


^ See for details Ref. [T3l and the review [T3l 
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Cdglap — 1 + 


0-s{Q^) ^( 1 ) 


IDGLAP = 



(26) 


Therefore, we can borrow from the DGLAP formulae the contributions which are missing in Eqs. 15I6II by adding 
Cdglap and 'Jdglap to the coefficient functions and anomalous dimensions of Eqs. m . Of course, both Cdglap 
and Jdglap contain also terms already included in Eqs. 15I6II . so in order to avoid a double counting these terms 
should be extracted from C^g, H^g. 

In order to do so, let us consider the region of a: ^ 1 where the effective values of lo in Eqs. 00 are large. In this 
region we can expand H^g and C^g into a series in l/co. Retaining the first two terms in each series and expressing 
them through A and we arrive at C^g = cj^ + O(a^), H^g = H^g + 0(af) with 


cSS’+ +1/2^1 . (22) 

~l±, ^ ^ ^ “(1/,^) [2/„ + 1] 2 + o(±)(„)[!/,., + i/2| , 

Now let us define the new coefficient functions C^g and new anomalous dimensions as follows: 


_ 

^NS — 


W - H 


(±) 

NS 



A(±) 

^NS 


^(±) 

'-'NS 


_ ^(±) 
'-'NS 


+ 1 + 


27r 


(28) 


These new, ’’implemented” coefficient functions and anomalous dimensions of Eq. include both the total 
resummation of the leading contributions and the DGLAP expressions in which as{Q'^) is replaced by A{uj). As shown 
in Ref. [^, this parametrization differs from the DGLAP one, as = as(fcl), at small x, though both parameterizations 
coincide when x is large. The main point is that the factorization of the phase space into transverse and longitudinal 
spaces used in DGLAP is a good approximation for large x only. The implemented coefficient functions C^^g and 
anomalous dimensions H^^g of Eq. 12811 can be used either in the standard DGLAP way by letting w be integer, i.e. 
uj = n, with n = 1, 2,... and thereby obtaining the moments of the non-singlet structure functions, or alternatively by 
using C^j^g,H^j^g in Eqs. 15I6II where uj is complex. In the latter case, Cdglap and Jdglap should first be continued 
analytically from the integer values ui = n. We have verified that the new formulae in Eq. 12811 provide a much better 
agreement with NLO DGLAP in the region of large x. 

We notice that our approach is qualitatively close to the one suggested in Ref. 0, however they are not quite 
identical: indeed they first differ in the treatment of Ug which is running in our picture in every Feynman graph 
involved, whereas in Ref. M the DL contributions are obtained from Refs, where as was kept fixed at an unknown 
scale. To make it running, in Ref. M the parametrization Ug = as{Q^) is suggested in the final formulae, which is 
of course incorrect at small x (see Refs. 0 for details). In addition, the single- logarithmic contributions were not 
accounted in Ref. M at all. 


VI. DISCUSSION AND CONCLUSIONS 

As is well known, the conventional DGLAP approach was originally suggested for the region of rather large values 
of X, where logarithms of gave the most important contributions. At the same time, DGLAP neglects the total 
resummation of logs of x. In our approach we have been able to account for those contributions to the non-singlet 
structure functions at small x. The results are presented in Eqs. 00. Similarly to DGLAP, the new anomalous 
dimensions 0 control the -evolution while new coefficient functions 0 evolve the initial quark densities from 
their empirical values at x ~ 1 to a: <C 1. 

The extrapolation of Eqs. 15I6II into the asymptotic region x ^ 0 leads to the Regge asymptotic formulae of 
Eqs. The extrapolation of the DGLAP expressions with the standard assumption of Eq. © for the initial 

quark distribution 5q also leads to the Regge behavior Eq. CT . However, our approach provides the Regge behavior 
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because of the total resummation of the leading logarithms of x whereas the DGLAP Regge asymptotics is generated 
by the phenomenological factor a;““ in the fit of Eq. © for 5q. This allows us to conclude that the singular part of 
the DGLAP input for the initial quark densities mimics the impact of total resummation and can be dropped when 
the resummation is done, simplifying the structure of the input densities to expressions regular in x. It also turns 
out that the DGLAP small-x asymptotics, using the fit of Eq. © for 5q for the non-sin glet , is more singular than 
our predictions. This contradicts the various analyses of the experimental data of Refs, [l^, which agree with our 
predictions. Our asymptotics also differ from DGLAP as far as the -dependence is concerned: indeed we predict 
an asymptotic scaling, where the non-singlet structure functions asymptotically depend on one argument ^ = Q'^ jx^ 
instead of two arguments and x, separately. 

Giving up the singular factor from the DGLAP fit Eq. © and using its regular part only, we have compared 
numerically our formulae Eqs. with the NLO DGLAP results. Such a comparison confirms that in absence of 
the singular terms the DGLAP small ^-dependence and Q^-dependence of the non-singlets are slower than ours. The 
results have been displayed in Figs. 1,2. 

Finally, in order to obtain an overall approach which is consistent and accurate in both regions of small and large 
X, we have suggested to combine our results and those of DGLAP by an appropriate definition of the anomalous 
dimensions and the coefficient functions. Indeed, the expressions of Eqs. 01 for the improved anomalous dimensions 
and coefficient functions account for both NLO DGLAP terms and our total resummation of the most important logs 
of X. Therefore, Eqs. can be used both at large and small x. Then the new inputs for 5q in this approach could 
be much simpler functions of x than in the standard DGLAP approach. 
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VIII. FIGURE CAPTIONS 

Fig. 1: The x -dependence for: Ri of Eq. ll^ (curve 1), i ?2 of Eq. (EHll (curve 2) and i ?3 of Eq. (1^ (curve 3). All 
curves correspond to = 10 GeV^. 

Fig. 2: The -dependence for R{Q'^) of Eq. |21. Gurves 1,2,3 correspond to a: = 10 10 10 respectively . 
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